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Abstract 

The four-quark condensates appearing in QCD sum rule are calculated nonperturbatively in a realistic 
QCD motivated extended Nambu-Jona-Lasinio model. The calculation is in the framework of the 
effective potential for local composite operators up to next-to-the-leading order in the 1/N expansion 
in which the non-factorized parts of the condensates are included. We show in this paper the possibility 
of explaining the phenomenological enhancement factor k = 2.4 , needed for fitting the data in 
the factorization approximation, as the contribution of the non-factorized parts of the four- quark 
condensates. 
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One of the interesting but difficult problems in nonperturbative QCD is the calculation of the quark 
and gluon condensates. The two-quark condensate breaks the chiral symmetry which explains the 
lightness of the low lying CP mesons, while all condensates are relevant to the well-known QCD sum rule 
[0. QCD sum rule incorporates the quark and gluon condensates with perturbative QCD calculations. 
For example, in the study of the sum rule for mesons, the vacuum expectation value (VEV) of two 
current operators < Tj (x)j B (0) > is considered, and the operator product expansion (OPE) gives 

i f d 4 xe^ < Tj A (x)j B (0) >=J2 C* B (q) < O n > , (1) 

where C AB {q) 's are perturbatively calculable Wilson coefficients and < O n > 's are various 
nonperturbative vacuum condensates. In practical applications, on the RHS of (1), only C AB {q) 's 
are calculated from perturbative QCD while the condensates < O n > 's are left as free parameters 
determined by experimental inputs. To reduce the number of free parameters, people often make the 
simple factorization approximation to express the four-quark condensate in terms of the two-quark 
condensate, i.e. 

< ^f^Fi/j > « F r < ^ > 2 , (2) 

where the Lorentz structure of F is 1 for scalar (S), ^75 for pseudoscalar (P), 7^ for vector 
(V), 757^ for axial-vector (A), for tensor (T) (for color and flavor non-singlet, the color and 

flavor group generators X a /2 , t{ should also be included ), and Fr is a constant depending on the 
structure of F and can be easily calculated Q. However, applications of QCD sum rule to various 
hadronic processes with the approximation (2), e.g. the sum rules for baryons [||] p-meson []|] and 
pseudoscalar meson [Q], show that the theoretical results are not quite good and phenomenologically 
the Fr < ipi/j > 2 term should be enlarged by a factor k ~ 3.6 for fitting the data @@[@]. Further 
study of the perturbative QCD corrections to the Wilson coefficient of the operator (-ipFiJ)) 2 shows 
that the 0(a s ) correction increases the Wilson coefficient by a factor of 1.5 |J which means that 
the enhancement factor should actually be k 2.4 . 

So far there is no theoretical explanation of the enhancement factor k 2.4 . A possible and 
natural conjecture is that the need for an enhancement factor might be due to the neglect of the non- 
factorized part of the four-quark condensate < tpFtjj'ipFip >nf in (2), but this is difficult to prove 
from the first principles of QCD. Several attempts have been made to discuss this problem. A careful 
calculation of the anomalous dimensions of the four-quark operators shows that the deviation from (2) 
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varies significantly with the renormalization scale, i.e. even if (2) is good at low energy it may fail at 
some high energy scale [||. In Ref.|7|, a general formalism of the effective potential for local composite 
operators in the framework of 1/N expansion was developed and was applied to calculate the four- 
fermion condensate < Tptptpip > nonperturbatively in the Gross-Neveu model. It is shown that (2) 
holds only in the N — > oo limit, i.e. < ipif)t/jifj >nf is nonvanishing to order-l/N . Up to order- 
1/N, for N = 3 , the numerical result shows that the ratio R =< 'tpipTpip >nf /[Fs < V'V' > 2 ] runs 
rapidly with the renormalization scale, and the minimal value of R is about 1/3 . This means that 
(2) is not a good approximation even at low energy. However, the Gross-Neveu model is a toy model 
which can not predict the value of k to compare with the experiment. In recent years, there have been 
several papers studying QCD motivated Nambu-Jona-Lasinio (NJL) models as effective theories for low 
lying hadrons^]. These models lead to quite successful phenomenological predictions, especially the 
extended Nambu-Jona-Lasinio (ENJL) model in Ref. Q, so that they may reflect some main features 
of nonperturbative QCD. Since the NJL model is a four-fermion interaction theory similar to the Gross- 
Neveu model, the method developed in Ref.[j^] can be directly applied to it. In this paper, we take 
the EJNL model of Ref.[|| and calculate the four-quark condensates appearing in the QCD sum rule 
nonperturbatively by using the method of Ref.|^]. The four-quark condensates are generally related to 

the sum rule in the form Cg < 'tpT^'^T^ > (cf. (1)). Once the condensates < tptp > and 

r 

< ipT^^Ttp > are calculated, the enhancement factor k can be obtained from 

J2 cl < > ]T cl < WW > NF 

- 1 + r ^ r ^ . (3) 



r r 

Here we sketch the main points of the method developed in Ref. 0], To calculate the two-quark and 
four-quark condensates, we consider the following generating functional 

Z[K 2 ,K 4 ] = expiW[K 2 , K 4 ] = i J VibV^expi J d 4 x{£ + ^ K\ W + ^ K\ W W 
+ P(K 2 ,K A )} , 

(4) 

where K\ and K\ are local external sources, and P{K 2 ,Ki) is a pure source term needed for 
satisfying the consistency condition requiring that no condensate effect occurs at the classical level |9|. 
Let W c be the classical value of W . The consistency condition is |§ 



Here, for simplicity, we have ignored the external sources for the single fields ip and tp in (4) since 
we are not calculating < ip > and < ip > which actually vanish. Define the classical fields £p 
and Hr , 

5K? = Sr ' 61$ = Sr + " r ' (6) 
and Hr + (1 — -^r)S| give < -ipip > and < ipTipipTip >nf , respectively, in the vacuum 
state. Once W[K2, K4] is obtained, we can make the Legendre transformation 

r[E,H] = W[K 2 ,K 4 ] - f d 4 x{Y,K T 2 Y, T + £ K\ ( + H r ) } , (7) 

r r 

to get the effective action r[S, E] and so the effective potential V e t t (S, 5) . The rules for calculating 



W[K 2 , K4] are as follows. Define the well-known propagator [l(|] 

5 2 S 



iD 1 (x,y;K 2 



5ip(y)5ip(x) 



(8) 

ip,i>=o 



where S is the action S = J (fx [C + £ K^ipTip + £ if[ ipTipipTip + P(K 2; if 4)] ■ £> 



r 

related to the physicsl propagator G by 

G -1 = v- 1 - n , (9) 

where II is the proper self-energy of the ip field. For models like the ENJL model, we can always 
consider a special proper self-energy diagram H s , composed of physical propagators and bare vertices, 
which is momentum-independent and is of the following structure [[/]] 

n s = -^6 r rA r , (10) 

r 

where 

5Wt 

Ars -«f = rr(rG) - 

and br is a constant determined by the model of interaction. In (11) Wl means the loop contribution 
to W . Define a new propagator G s , 

G- 1 = V- 1 - U s . (12) 

In terms of these notations, the formula for W[K 2 , Ka\ derived in Ref. 0] is written as 

W[K 2 ,K 4 ] = I d A x{C + J2 K 2^ + X)- K 4#W I > + P(K 2 ,K 4 )} - iTrlniG' 1 

r r r f . , (13) 

— / d 4 x ^6rAp — i < 0\expi / d 4 x£j\0 >'p2Pim s ) > 



where Cj is the interaction Lagrangian, and the last term means the sum of all the partially two-particle 
irreducible diagrams with respect to II S (with propagators G s 's and vertices in Cj ) defined in Ref. 
[[?]] except the type of diagrams shown in Fig. (4c) in Ref. 07||. The first term in (13) is the tree level 
contribution, while other terms are contributed by loop diagrams ( Wl )■ In (13), A serves as an 
unknown constant which can be further solved from (11). Then the effective potential V e ff can be 
obtained from (13) and (7). The condensates can then be determined by the minimum of V e ff ■ 

Actually, once W is obtained, we can also calculate the condensates simply from (6) by taking 
Ki = K4 = without knowing V e ff . The solutions of (6) are the same as those obtained from the 
extremum condition of Veff as it should be[^]. To pick up the true condensates (minimum of V e ff ), 
one needs only calculate the values of V e ff f° r these solutions. This is much easier than calculating 
the complete function V e jy(£,E) for arbitrary £ and S. 

The Lagrangian in the ENJL model in Ref. []|] is 

£QCD = Cq X CD + CffjL + £ V N j L , (14) 

where A x is a momentum cut-off, I^qcd ' s tne QCD Lagrangian for momentum below A x , and 

^ = ^E(«£)(to) 

c X ab 

2n 2 G s r , 7n2 



[(^V) 2 + Wi75^) 2 ] (15) 



2n 2Q s N h l 

+ AT AT A 2 E [OMiVO(V^) + (^75^X^75^)] ) 

£ v n a jl = -^^E^wixfev^) + ^w b R )^Wr R )] 

c X ab 

= ~ N N A 2 KV^VX^VO + ($7&W)($Tb7/^0] (16) 
4ir 2 G N2f ~ l 

~ AT AT I2 E [^l^M^lM) + {^l^)(^lM)\ ■ 
N c N f A X i=l 

in which Gs , Gy are two coupling constants treated as free parameters, Gs = NfGs , Gy = NfGy , 
and the flavor group generator t{ is normalized as ti = Xi/\/2 , i = 1, 2, • • • , Nj — 1 . It is argued in 
Ref. [U that the NJL type Lagrangian Cff P JL + C^ul mav De understood as coming from integrating 
out the high momentum modes of quarks and gluons in the fundamental theory of QCD and is regarded 
as the main part in Cqcd that is responsible for the chiral symmetry breaking, while ^qcd ' s 
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responsible for the low-energy gluonic corrections to the broken chiral symmetry state. This supports 
the idea of the chiral quark model [11]. Even without the gluonic corrections from £qcd , the NJL 
type Lagrangian can give rather good phenomenological predictions for low energy hadron physics Q. 
Since we are not aiming at accurate calculations, we simply neglect the gluonic corrections in this paper 
and keep only Cf/jL + ^tvjl as tne effective interactions between quarks, i.e. the Lagrangian C in 
(4) is taken to be 

C = $ipil> + Cf f P jL + t%jL- (17) 
The coefficients br in (10) and (13) are now 



b S = b P = AT AT , 2 > b S = b p- 



8tt 2 < 

V — U A 



b i = b i _ _ 8tt 2 G v hU Lt . 8tt 2 G v 



(18) 



and the pure source term P{K2,K^) obtained from (5) is 

P(K 2 , KA = -y( ^ -) . (19) 

There are three independent free parameters Gs , Gy , and A x in the ENJL model in Ref. ||. 
They are related to the constituent chiral quark mass Mq , the quark axial-vector coupling constant 
gA , and the vector meson mass My by the following relations^] 

l/G s = (M Q /A x ) 2 r(-1, (M Q /A X ) 2 )(1 + 7 _i) , (20) 

9A = 1 + 4G v (M Q /A x mo, (M Q /A X )2)(1 + 701 ) ' (21) 

A 2 = 1m v G v T(0, (Mq/A x ) 2 )(1 + 703 ) , (22) 

where T{n — 2, (Mq / A x ) 2 ) is the incomplete gamma function ||, and 7 _i , 701 , and 703 are 
gluonic corrections which are actually not large and have been neglected in our Lagrangian (17). Apart 
from the gluonic corrections, the rest part of the formulae used to fit the data in Ref. [Q] are all of the 
leading order in 1/N C expansion. Actually Mq , A x , and are taken as the input parameters 
to fit the data in Ref. ||. There are five different ways of fitting the data presented in Ref. |8|], which 
determine different sets of values of Mq , A x , and gA , and the predictions are all successful. The 



simplest one of the fits is their Fit 4 in which g^ = 1 . From (21) we see that this corresponds to 
Gy = which will simplify our calculation. Moreover, it is explained by Weinberg [^] that should 
actually be unity in the constituent quark model. In view of this, we take in this paper the set Fit 4 
in Ref. [§, and by means of (20)-(22), it leads to 

G s = 1.19, G v = 0, A x = 667 MeV . (23) 

In the present model, N c = Nj = 3 . Let N = N c = Nf . Our dynamical calculation of the 
condensates is in the l/N expasion. In doing so, we should take a proper limit of N —> oo such that 
the three terms in (17) are of the same order in this limit, i.e. both the propagation of the quark field 
and all the NJL type interactions are included. This means that (cf. (15)-(17)) we should treat the 
flavor singlet coupling constants Gs , Gy and the flavor non-singlet coupling constants Gs , Gy 
as independent finite parameters as N — > oo . After obtaining the condensates to certain order in 
l/N expansion, we then take numerically N = 3 and Gs = 3Gg , Gy = 3Gy when calculating 
the enhancement factor k . In order to make comparison with the formulae in Ref. ||, we define the 
constituent chiral quark mass as 

Mq = - b s < W > , (24) 

which corresponds to the mean-field approximation interpretation of Mq. Using the method described 
above, similar to the calculation presented in Ref. 0], we obtain the following results. 

In the large- iV limit, W[K2, K&] is contributed by chain diagrams with open ends. Then the gap 
equation (eq.(6) with K2 = K4 = ) to this order reads 

<*>= _, tr ^_^_ = __2r ( -i,(W), (25) 

< ^75^ > = < iplfiip > = < ^757^ > = . 
In (25) the f A means that the A x -regulator used in Ref. H is taken for the momentum integration. 
In the chiral limit < xptp >= N < uu > , thus we have 

< uu > = - -?L Mq r(-l, (M Q /A X ) 2 ) . (26) 

which coincides with the result in Ref. ||. It is also easy to see that, in the large- N limit, [SW /SK^]k 2 ,k 4 
Sp , so that [^r]K 2 ,K 4 =o vanishes. Since in this limit Fp = 1 , \^v\K 2 ,K i =Q is just the non- factorized 
part of the four-quark condensate. Thus we have 

< vry^rv> >nf= o (27) 



in the large-iV limit. The intuitive reason for the factorization of the four-quark condensate in the 
N — > oo limit is as follows. Diagrammatically, the four-quark condensate is obtained by removing a 
node from the chain diagram for W , which releases four quark lines. For open-end chains, the removal 
of a node causes two disconnected chains , each with two quark lines, and this is just the factorization 
of the four-quark condensate. 

To order- 1/N , we should further take into account the P2PI(n s ) diagrams shown in Fig. 6 in Ref. 
|0| , which include the closed chain diagrams. Removing a node from a closed chain does not cause two 
disconnected chains, so that factorization breaks down to this order. Up to Next-to-the- leading order 
in 1/N expansion, we see from (3) that only an order- 1/N calculation of < tpTipipFtp >nf is 
further needed. 

In practical applications, the currents in (1) are usually color octets. Similar to Ref.[[/]], by applying 
the method described above, we obtain the following up to order-l/iV results of the non-factorized 
parts of various relevant four-quark condensates. 

< uX a uuXaU >NF = {F S /N)(S NF - P NF - A NF ) + 0{1/N 2 ) , 

< uX a ddX a u > NF = (F S /N)S NF + 0{1/N 2 ) , 

< uX a udX a d >nf= + 0(l/iV 2 ) , 



< ui-f 5 X a uui"f 5 X a u > NF = {F P /N)(S NF - P NF + A NF ) + 0(1/N 2 ) 

< uij 5 X a ddij 5 X a u >nf= {F P /N)S NF + 0(1/N 2 ) , 

< uij5X a udij5X a d >nf= + 0(1/N 2 ) , 



< u^\ a uu^\ a u > NF = {F V /N)(S NF + P NF - -A NF ) + 0(1/N 2 ), 

< u^X a dd^X a u > NF = (F V /N)S NF + 0(1/N 2 ) , 

< U^XaUdju^ad >NF = + 0(1/ N 2 ) , 



(29) 



(30) 



< uj 5 ^X a uu-f 5 ^X a u > NF = (F A /N){S NF + P NF + -A NF ) + 0(1/N 2 ), 

< U'J^Xadd^J^XaU >NF= ~^ S NF + 0(1/N 2 ) , 

< u^X a ud^^X a d > NF = + 0(1/N 2 ) , 
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< ua^Kuua^Ku > NF = (F T /N)(S NF - P NF ) + 0(1/N 2 ) , 

< uo^\ a ddo^\ a u >nf= (F T /N)S NF + 0(l/iV 2 ) , 

< ua^ u X a uda^ v \ a d > NF = + 0(l/iV 2 ) , 



where 



and 



F s = -F P = -1 + 0(l/N 2 ), 

Fy = — Fa = —2 + 0(l/iV 2 ) , (33) 
F T = -6 +0(l/iV 2 ), 

ivr 2 G 5 /■ d 4 g T 2 (q 2 ,M Q ) 



Snf 



4 ^ A x (2*)' i _ ^|Tp(^, Mq) ' (34) 



4iV A 2 y Ax (2vt)4 1 _^ T5((? 2 )Mq) 



i^ 2 G 5 /■ A T P (q 2 ,M Q ) 



A 



NF 



iir 2 G s f d*q q- 2 Tp A (q 2 ,M Q )T A p(q 2 ,M Q ) 



in which 



Ts(q 2 ,M Q )=^Nl m 



4 iVA| 7 Ax (2vr)4 1-^T p (^Mq) 



d 4 p . 1 



A x (2vr)4 L ^ + ( / l /2-M Q {>- i/2 - Mq 1 
T(-l, (-« ) 2 ) - (2g 2 - 8M 2 ) I^ 2 , M Q ) , 



T P (q 2 ,M Q ) = AiN f 
J A 



d A p , 1 1 



A x (2vr)4 L '*t + 4/2-M Q '^-4/2-Mq 1 
VM ir(-l,(^) 2 )- V/(, 2 ,Mq), 05) 



- 2 ^ *' V A X 

dp 



Tp,(, 2 ,M ) = -W,Mq) = ^L^[i, 5 l5 4 



A x (2^) 4 L '°t + 4/2-M Q '^-^/2-Mq j 

= -iMg^ 2 /^ 2 , Mq), 

iV f 1 q 2 (x 2 -x)+M'k 

In (35) tr is the trace in the Dirac spinor space. With all these, we can calculate the enhancement 
factor k for practical processes. 

Take the p-meson sum rule as an example. The relevant four-quark condensates are [Q 

C 6 < 6 > = 6vr 3 a s {< {u^lbKu - d^-f 5 X a d) 2 > 

+ ^ < {u^X a U + Jy*Aad) ^ V'T/.Aa^ >}• ^ 

u,d,s 
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Substituting (36) into (3), and using the results in (30) and (31), we obtain the enhancement factor 
for the /9-meson sum rule 

1 SnF + PNF . 11 Anf r^, 1 \ /o-r-v 

K = 1 + 3 <uu>* + 42 7^ + ° ( A^ " (37) 

When doing numerical calculations, we should bear in mind that the numbers in (23) are obtained 
from fitting the data with the large-iV c formulae together with gluonic corrections, while our present 
calculation is up to order-l/iV without gluonic corrections. Thus we should not take the numbers in 
(23) so seriously when calculating the value of k . From (36) we see that the values of Snf, Pnf 
and Anf are only related to the two parameters Gs and Mq/A^ which are further related to 
each other by (20). Therefore our k is a function of only one parameter Gs ■ In Table 1, we list 
the values of Snf, Pnf and Anf for several values of Gs around the number given in (23). The 
corresponding numbers of k in the p-meson sum rule are listed in Table 2. We see that in the small 
Gs region k is a sensitive function of Gs , and k = 2.4 is included in this region. At present 
Gs is not determined accuratly, so that we still cannot predict quantitatively if k is precisely 2.4. 
However, due to this sensitivity, the possibility of having k = 2.4 mainly depends on the positivity 
of the order-l/iV contribution in (37) (the sum of the second and third terms in (37)) which is not 
obvious in general. Our calculation shows that it is really the case. 

Inview of the fact that k is sensitive to the value of Gs which is not yet well determined, our 
conclusion in this paper is that the explanation of k = 2.4 as the contributions of the non-factorized 
parts of the four-quark condensates is possible. 
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Table 1 Values of Snf , Pnf , and Anf (in the unit of <uu> 2 ) for given values of 



G s 


2.0 


1.5 


1.3 


1.25 


1.2 


1.15 


1.1 


Snf 


0.00069 


0.0028 


0.0079 


0.011 


0.016 


0.026 


0.047 


Pnf 


0.010 


0.017 


0.026 


0.030 


0.037 


0.048 


0.072 



A NF -0.00054 -0.00077 -0.00094 -0.00099 -0.0011 -0.0011 -0.0012 



G s 


1.08 


1.05 


1.05 


1.02 


1.01 


1.009 


1.008 


Snf 


0.065 


0.12 


0.23 


0.37 


0.85 


0.96 


1.1 


Pnf 


0.090 


0.15 


0.26 


0.40 


0.88 


0.99 


1.1 


Anf 


-0.0012 


-0.0013 


-0.0013 


-0.0014 


-0.0014 


-0.0014 


-0.0014 




G s 


1.007 


1.006 


1.005 


1.004 


1.003 


1.002 


1.001 


Snf 


1.3 


1.5 


1.9 


2.4 


3.4 


5.3 


12 


Pnf 


1.3 


1.6 


1.9 


2.5 


3.4 


5.3 


12 


Anf 


-0.0014 


-0.0014 


-0.0014 


-0.0014 


-0.0014 


-0.0014 


-0.0014 



Table 2 Values of the enhancement factor k ( eq.(37)) for given values of Gs ■ 



G s 


2.0 


1.5 


1.3 


1.25 


1.2 


1.15 


1.1 




1.0035 


1.0065 


1.011 


1.014 


1.017 


1.024 


1.039 




G s 


1.08 


1.05 


1.03 


1.02 


1.01 


1.009 


1.008 


fv 


1.051 


1.089 


1.16 


1.26 


1.57 


1.65 


1.74 




G s 


1.007 


1.006 


1.005 


1.004 


1.003 


1.002 


1.001 


K 


1.86 


2.0 


2.3 


2.6 


3.3 


4.6 


8.7 
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